Waiting time analysis of foreign currency exchange 
rates: Beyond the renewal-reward theorem 
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Abstract — We evaluate the average waiting time between 
observing the price of financial marliets and the next price 
change, especially in an on-line foreign exchange trading service 
for individual customers via the internet. Basic technical idea 
of our present work is dependent on the so-called renewal- 
reward theorem. Assuming that stochastic processes of the 
market price changes could be regarded as a renewal process, 
we use the theorem to calculate the average waiting time of 
the process. In the conventional derivation of the theorem, it is 
apparently hard to evaluate the higher order moments of the 
waiting time. To overcome this type of difficulties, we attempt to 
derive the waiting time distribution fl{s) directly for arbitrary 
time interval distribution (first passage time distribution) of the 
stochastic process Pw{t) and observation time distribution Po{t) 
of customers. Our analysis enables us to evaluate not only the 
first moment (the average waiting time) but also any order of the 
higher moments of the waiting time. Moreover, in our formalism, 
it is possible to model the observation of the price on the internet 
by the customers in terms of the observation time distribution 
Po{t). We apply our analysis to the stochastic process of the 
on-line foreign exchange rate for individual customers from the 
Sony bank and compare the moments with the empirical data 
analysis. 

I. Introduction 

Fluctuation has an important role in lots of phenomena 
appearing in our real world. For instance, a kind of magnetic 
alloy possesses magnetism in low temperature, whereas in 
high temperature, it losses the magnetism due to thermal 
fluctuation acting on each spin (a tiny magnet in atomic scale 
length). Thus, the large system undergoes a phase transition 
at some critical temperature and the transition occurs due 
to cooperative behavior of huge number of spins, to put it 
other way, due to competition between thermal fluctuation and 
exchange interaction between spin pairs making them point to 
the same direction [1]. 

To understand these kinds of macroscopic properties or 
collective behavior of the system from the microscopic point 
of view, statistical mechanics provides a good tool. Actually, 
statistical mechanics has been applied to various research 
subjects in which fluctuation is an essential key point, such 
as information processing [2] or economics including finan- 
cial markets and game theory [3]. Especially, application of 
statistical-mechanical tools to economics, data analysis of 
financial markets — what we call econophysics — is one of 
the most developing research fields [4], [5], [6]. Financial data 



have attracted a lot of attentions of physicists as informative 
materials to investigate the macroscopic behavior of the mar- 
kets from the microscopic statistical properties [4], [5], [6]. 
Some of these studies are restricted to the stochastic variables 
of the price changes (returns) and most of them is specified 
by a key word, that is to say, /of tails of the distributions [4]. 
However, the distribution of time intervals also might have 
important information about the markets and it is worth while 
for us to investigate these properties extensively [7], [8], [9], 
[10], [11], [12]. 

Such kinds of fluctuation in time intervals between events 
are not special phenomena in price changes in financial 
markets but very common in science. In fact, it is well-known 
that spike train of a single neuron in real brain is time series 
in which the difference between successive two spikes is not 
constant but fluctuated. This stochastic process specified by the 
so-called Inter-Spike Intervals (ISI) is one of such examples 
[13], [14]. The average of the ISI is about a few milli-second 
and the distribution of the intervals is well-described by the 
Gamma distribution [14]. 

On the other hand, in financial markets, for instance, 
the time intervals of two consecutive transactions of BUND 
futures (BUND is the German word for bond) and BTP 
futures (BTP is the middle and long term Italian Government 
bonds with fixed interest rates) traded at LIFFE (LIFFE 
stands for London International Financial Futures and Options 
Exchange) are ^ 10 seconds and are well-fitted by the so- 
called Mittag-Leffler function [8], [9], [10]. The Mittag-Leffler 
function behaves as a stretched exponential distribution for 
short time interval regime, whereas for the long time interval 
regime, the function has a power-law tails. Thus, the behavior 
of the distribution described by the Mittag-Leffler function is 
changed from the stretched exponential to the power-law at 
some critical point [15]. However, it is non-trivial to confirm 
if the Mittag-Leffler function supports any other kind of market 
data, for example, the market data filtered by some rate 
window. 

As such market data, the Sony bank USD/JPY exchange rate 
[16], which is the rate for individual customers of the Sony 
bank in their on-line foreign exchange trading service, is a 
good example to be checked by the Mittag-Leffler function. 
Actually, our preliminary results imply that the Mittag-Leffler 



TABLE I 

Typical three examples with fluctuation between the events. 





ISI 


BUND future 


Sony bank rate 


Average time interval 


~ 3 [ms] 


~ 10 [s] 


~ 20 [min] 


PDF 


Gamma 


Mittag-Leffler 


WeibuU 



function does not support the Sony bank rate [17]. The Sony 
bank rate has ~ 20 minutes [18] as the average time interval 
which is extremely longer than the other market rate as the 
BUND future. This is because the Sony back rate can be 
regarded as the so-called first passage process [19], [20], 
[21], [22], [23], [24] of the raw market data. In Table U we 
list the average time intervals and the probability distribution 
function (PDF) that describes the data with fluctuation between 
the events for typical three examples, namely, the ISl, the 
BUND future and the Sony bank rate. From this table, an 
important question might be arisen. Namely, how long do the 
customers of the Sony bank should wait between observing 
the price and the next price change? This type of question 
is never occurred in the case of the ISl or the BUND future 
because the average time intervals are too short to evaluate 
such informative measure. 

Obviously, for the customers, an important (relevant) quan- 
tity is sometimes a waiting time rather than a time interval 
between the rate changes. The waiting time we mentioned here 
is defined by the time for the customers to wait until the next 
price change since they try to observe it on the World Wide 
Web for example [16]. If the sequence of the time intervals 
has some correlations and the customers observe the rate at 
random on the time axis, the distribution of the waiting time 
is no longer identical to the distribution of the time intervals. 
In the previous studies [25], [26], [27], we assumed that the 
time intervals of the Sony bank USD/JPY exchange rates 
might follow a WeibuU distribution and evaluated the average 
waiting time by means of the renewal -reward theorem [28], 
[29]. However, the conventional renewal -reward theorem is 
restricted to the case in which the customers observe the rate 
at random on the time axis and it is hard to extend the theorem 
to the situation in which the time for the customers to observe 
the rates obeys some arbitrary distributions. 

To make these problems and difficulties clear, in this paper, 
we introduce a different way from the renewal-reward theorem 
to evaluate the higher-order moments of the waiting time 
for arbitrary time interval distribution of the price changes 
and observation time distribution by directly deriving the 
waiting time distribution. We first show that the result of 
the renewal-reward theorem [27] is recovered from our new 
formalism. Then, it becomes clear that our formulation is more 
general than the conventional renewal-reward theorem. As an 
advantage of our approach over the renewal-reward theorem, 
we can evaluate the higher-order moments of the waiting 
time, and moreover, it becomes possible to consider various 
situations in which the customers observe the rate according 
to arbitrary distribution of time. 

This paper is organized as follows. In the next section HIl 



we introduce the Sony bank rate [16] which is generated from 
the high-frequency foreign exchange market rate via the rate 
window with width 2e yen (e = 0.1 yen for the Sony bank). 
In section |llll we summarize a series of our previous studies 
related to the present paper To understand the mechanism of 
the Sony bank rates as the first passage process [19], [20] of 
the raw market rate, in section IIVI we carry out the computer 
simulations by making use of the GARCH (Generalized Au- 
toRegressive Conditional Heteroscedasicity) model [30], [31], 
[32] with discrete returns as the time series behind the Sony 
bank rates. The effect of the rate window on the currency 
exchange rates is revealed. In the next section [V] we explain 
our method to derive the waiting time distributions. We show 
that our treatment reproduces the result by the renewal-reward 
theorem. We also evaluate the deviation around the average 
waiting time for the Weibull first passage time distribution 
and uniform observation time distribution. We find that the 
resultant standard deviation is the same order as the average 
waiting time. We test our analysis for several cases of the 
observation time distributions and calculate the higher-order 
moments. The last section [Vl] is concluding remarks. 

11. The Sony bank rates as a first passage process 

The Sony bank rate (see Fig. [TJ we shall deal with in this 
paper is the rate for individual customers of the Sony bank 
[16] in their on-line foreign exchange trading service via the 
internet. If the USD/JPY market rate changes by greater or 

Sony bank rate 




ticks 



Fig. 1. Behavior of the Sony bank rates. 

equal to 0.1 yen, the Sony bank USD/JPY exchange rate is 
updated to the market rate. In this sense, the Sony bank rate 
can be regarded as a kind of first passage processes [19], 
[20], [21], [22], [23], [24]. In Fig.|2l we show the mechanism 
of generating the Sony bank rate from the market rate (This 
process is sometimes refereed to as first exit process [33]). 
As shown in the figure, the difference between successive 
two points in the Sony bank rate becomes longer than the 
time intervals of the market rates. In Table [III we show 
several data concerning the Sony bank USD/JPY rate vs. 
tick-by-tick data by Bloomberg for USD/JPY rate. It is non- 
trivial problem to ask what kind of distribution is suitable 




market rate 

• Sony bank rate 



first passage time 



Fig. 2. An illustration of generating the tiltered rate (black circle) by the 
rate window with width 2e (shaded area) from the market rate (solid line). 

TABLE II 

The Sony bank USD/JPY exchange rate vs. tick-by- tick data 
for usd/jpy exchange rate. 





Sony bank rate 


tick-by-tick data 


# of data a day 


~ 70 


- 10,000 


The smallest price change 


0.1 yen 


0.01 yen 


Average interval between data 


~ 20 minutes 


~ 7 seconds 



to explain the distribution of the first passage time. For this 
problem, we attempted to check several statistics from both 
analytical and empirical points of view under the assumption 
that the first passage time might obey a non-exponential 
Weibull distribution [25], [26], [27]. We found that the data is 
well explained by a Weibull distribution. This fact means that 
the difference between successive Sony bank rate changes is 
fluctuated and has some memories. 

III. Previous related results 

Before we start to explain what we attempt to do in the 
present paper, we shortly summarize the related results which 
were akeady obtained by the present authors. 

• Sazuka (2006): We checked that the Sony bank rate is 
well-described by a Weibull distribution by some empir- 
ical data analysis (Weibull paper analysis, evaluation of 
Kullback-Leibuler, Hellinger divergence measures) [25], 
[26]. 

• Inoue and Sazuka (2006): We showed analytically that 
the crossover between non-Gaussian Levy regime to 
Gaussian regime is observed even in the first passage 
process (which is the same process of the Sony bank 
rates) for a truncated Levy flight (the so-called KoBoL 
process [34], [35], [36] in mathematics or mathematical 
finance) [37]. 

• Inoue and Sazuka (2006): We introduced queueing the- 
oretical approach into the analysis of the Sony bank rate 
and evaluated average waiting time including expected 
returns. We also carried out computer simulations by 



using the GARCH model to investigate the effect of the 
rate window of the Sony bank [27]. 

• Sazuka (2006): We observed a phase transition between 
a Weibull distribution to a power-law distribution at some 
critical time from the empirical data analysis of the Sony 
bank rate [26]. 

• Sazuka and Inoue (2007): We introduced the Gini index 
to evaluate to what extent a Weibull distribution is well- 
fitted to explain the behavior of the first passage process 
of the Sony bank rate. The analytical evaluation and the 
empirical data analysis gave quite similar results [38], 
[17]. 

Then, we focus on the following two points. 

> The effect of discreteness of returns in computer simula- 
tion by means of the GARCH model. 

> Generalization of the renewal-reward theorem to calculate 
the higher order moments of the average waiting time or 
to evaluate the average waiting time for the case in which 
the observation time distribution of the trader is explicitly 
given. 

This paper is intended as an investigation of these two points. 
IV. Computer simulations by the GARCH model 

WITH discrete returns 

In the previous studies [27], we carried out computer 
simulations of the Sony bank USD/JPY exchange rates by 
assuming that the raw market data might be well-described by 
the GARCH model [30], [31], [32] in which the successive 
time intervals At obey a Weibull distribution specified by a 
single parameter mp, namely. 



P{M) = 



niQ (At)™"- 



■ exp 



(Aty 



(1) 



In [27], we set a = 1 for a scaling parameter for simplicity to 
handle. We investigated the stochastic process generated by the 
rate window and estimated the distribution of time intervals of 
the stochastic process filtered by the rate window. Then, we 
assumed that the output of the filter having rate window with 
width 2e also follows a Weibull distribution with to(^ toq) 
and estimated m by means of Weibull paper analysis. The 
plot m-mp was obtained and the effect of the rate window on 
the market rate became clear. However, in those simulations, 
we assumed that the smallest price change (return) might take 
continuous values for simplicity to carry out the simulations. 
As we see in Table 1, the smallest price change unit of the tick- 
by-tick market data is 0.01 yen. Therefore, we should modify 
our computer simulations of the GARCH model by taking into 
account the discreteness of the price change. In this section, 
we show the result of the computer simulations. To this end, 
we deal with here the foUowing modified GARCH( 1,1) model 
[27]: 



= ao 



aiX. 



t-At 



At 



(2) 
(3) 



where the time interval At obeys a Weibull distribution with 
parameters (m, a), namely, 

m(Ai)"-i 



P{At) 



■ exp 



[At)- 



(4) 



and A/'(0, tTj) in this expression stands for a Gaussian with 
zero-mean and time dependent variance af. For simplicity, 
we set the parameter a as a = 1. We should notice that the 
GARCH(1,1) model has the variance cr^ after observing on a 
long time intervals t ^ oo and leads to 

2 "0 



1 — Qfl — /3l 

To make the return Xt discrete, for each time 
regenerate Xt by using the next map: 

Xt = *A(^t) = Aceil(A-iXt) 



(5) 



step, we 



pdf of Ft = Xt — Xt-At for the above discrete GARCH model 
with a parameter set: {ao, ai, (3i) = (0.4, 0.3, 0.3) which gives 
the kurtosis k = 4.17 and variance cr^ = 1 in the limit of 
< ^ oo, and A = 0.1 and e = a = 1. We should keep in 
mind that the ratio between the width of the rate window and 
the minimal change of the rate, namely, e/A ~ 10 is the same 
as that of the Sony bank rate. From this figure, we recognize 
that the return Y actually takes discrete values as expected. 

To investigate the effect of the rate window with width 2e, 
we estimate the Weibull parameter m of the output sequence 
from the rate window by means of the so-called Weibull 
paper analysis [25], [26], [27] for the cumulative Weibull 
distribution. From Fig. |4] we find that for the case of mo — 



(6) 



where the function ceil(a;) is defined as the smallest integer no 
less than x. The parameter A appearing in (|6|l means the length 
of the minimal variation of the return. In Fig. |3] we plot the 
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Fig. 4. The Weibull paper for the case of mo = 0.6 (upper panel) and 1.2 
(lower panel). We choose A = 0.1. 

0.6, 1.2, the first passage time distribution also obeys a Weibull 
distribution with a parameter m which is different from mg. 
In Fig. |5] we plot the relation between mo and m for several 
values of A. From this figure, we find that the relation for 
c- Tu jf f 1^ J i J- . . . the discrete cases with A = 0.1 and 0.2 are almost same as 

rig. 3. Ihe pat of the CjARLH(l,l) model with discrete returns. We set 

A = 0.1 (upper panel) and A = 0.2 (lower panel). The inset is behavior the relation for the continuous case A = reported in our 
around small Y regime. previous papers [27]. Meaningful differences are observed if 
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Fig. 5. The relation between mo and m for the case of A = 0,0.1,0.2 
and 0.5. For each A, we carried out 10 independent runs of simulations. 



we increase the value of A up to 0.5. This resuh provides 
us a justification of our previous GARCH modeHng [27] of 
the market rate to simulate the Sony bank USD/JPY exchange 
rates on the assumption that the minimum price change could 
take infinitesimal values. 

V. Derivation of the waiting time distribution 

In the previous studies [27], we evaluated the average 
waiting time for the customers to wait by the next price 
change since they attempt to observe the price by making 
use of the renewal-reward theorem [28], [29]. However, the 
theorem itself is obviously restricted to deriving only the first 
moment of the waiting time. From this reason, it is very 
hard to evaluate, for instance, the standard deviation from the 
average waiting time within the framework of the theorem 
(for example, see the proof of the theorem provided in [29]). 
Thus, we need another procedure to calculate it without the 
conventional derivation of the renewal-reward theorem. In this 
section, we directly derive the distribution of the waiting time. 
Our approach here enables us to evaluate not only the first 
moment of the waiting time but also any order of the moment. 
This section is a core part of this paper. 

A. The probability distribution of the waiting time 

We first derive the probability distribution function of the 
waiting time s. Then, let us suppose that the difference 
between successive two points of the Sony bank rate change, 
namely, the first passage time r follows Pw{t). Then, the 
customers observe the rate in time i (0 < t < r) that should be 
measured from the point at which the rate changes previously. 
In Fig. |6] we show the relation among these points r, t and s 
in time axis. The waiting time of the customers is naturally 
defined hy s = t — t. We should notice that the distribution 
Q,{s) is written in terms of the first passage time distribution 
Pvi/{T) and the observation time distribution Po(t) of the 



Fig. 6. The relation these points t, t and s in time axis. The first passage 
time r is given by r = t{i + 1) — t{i). The observation time is measured 
from the point t{i). 



customers as 

^{s) (X dT dtQ{s\T,t)Po{t)Pw{T). (7) 

Jo Jo 

Obviously, probability Q{s\T,t) that the waiting time takes s 
provided that the observation time and the first passage time 
are given as t and t, respectively, is given as 



Q{s\T,t) = 6{s-T + t) 



(8) 



with the delta function S{x). Taking into account the normal- 
ization constant of fl{s), we have 



n(s) 



drPwir) /; dt 6{s-T + t)Po{t) 
ds drPw (r ) /; dt5{s-T + t)Po {t) 



(9) 



where t denotes the observation time for the customers. We 
should notice that the result of the renewal-reward theorem : 
w = (s) — E{t'^)/2E{t) (see for example [29]) is recovered 
by inserting the uniform observation time distribution Pq (t) — 
1 into the above expression as 



(s) = / dssQ{s) 



dss drPwir) 
ds drPwiT) 



Io^fA-}d-CdrPwir) 
(1/2) J^°^s'Pwis)ds ^ E{r^) 
j;^sPw{s)ds 2E(t) 



(10) 



where we defined the rt-th moment of the first passage time 
E{t-^) by 



£;fr") 



/•OO 

/ dss'^Pwis). (11) 
Jo 



More generally, we may set Po{t)- For this general form of 
the observation time distribution, the probability distribution 
of the waiting time s is given as follows. 



n{s) 



dTPwiT)PoiT - S) 
ds dTP^{T)Po{T - 

dTPw{r)Po{T ~ s) 



E{t) - S, 



where we defined (5„ by 

dsi 







n 



Pw{t)- 



ds 



(12) 



(13) 



By using the same way as the derivation of the distribution 
ri(s), we easily obtained the first two moments of the waiting 
time distribution as 



is) 



E{t) - 5, 
and the standard deviation leads to 



E{t) - 5, 



'{AE{T^)E{T)-iE{T^)} + Gs, 



(52,(53 



12{E{t)-5^Y 



M = 



-ASiE{t^) 
1262E{t^) 



' 1263E{t) 
I261S3 - \2E\ 



where we defined 



dss"rj(s). 



(14) 



(15) 



(16) 



(17) 



Thus, this probability distribution r2(s) enables us to evaluate 
any order of the moments for the waiting time. 

In following, we evaluate the average waiting time and the 
deviation around the average for typical two cases, namely, 
Po(t) = 1 and Po{t) = c-*/^». 

1) The customers' observation follows Po{t) = I-' We 
first consider the case of Po{t) — 1. This case corresponds 
to the result obtained by the renewal-reward theorem [27]. 
Obviously, we find that (5„ = holds for arbitrary integer n. 
Thus, the waiting time distribution f2(s) leads to 



E{t) ■ 



(18) 



Then, the average waiting time and the deviation around the 
value lead to 



w 



E(r') ^_j AE{T^)E{TyiE{T^y ^ ^^^^ 



2E{Ty 



12E{tY 



For a Weibull distribution having the parameters to, a, the 
above results are rewritten by 

m = 4^^TTTTT (20) 



w 



r(-) 

ai/"Vr(l/TO)r(3/m) - r(2/TO)2 

^ ^ r(i/TO) 

where we defined the Gamma function as 



T{x) 



dte 



(21) 



(22) 



(23) 



It is important for us to notice that for an exponential distri- 
bution TO = 1, we have u; = a = a by taking into account 
the fact that T{n) = {n — 1)!. Moreover, the average waiting 
time w is identical to the average time interval E{t) since 
w = E{t^)/2E{t) = E{t) holds if and only if m = 1 (The 
rate changes follow a Poisson arrival process). These results 
are already obtained in our previous studies [27]. In Fig. |7] 



we plot the distribution D,{s) and the standard deviation a 
for a Weibull distribution. Especially, for the Sony bank case 
a = 49.6345, to = 0.585, we find a = 60.2284 minutes. 




Fig. 7. The distribution of waiting time for a Weibull distributing f2(s) with 
a = 1 and m = 0.59, 1 and 2 (upper panel). In the lower panel, we plot 
the standard deviation a for a Weibull distribution with a = 49.6345 as a 
function of m. 

2) The customers' observation follows Po{t) = e"*/"^".- We 
test the other observation time distribution of the customers. 
For instance, we might choose Po{t) — e"*/^" on the 
assumption that the customers might observe the rate more 
frequently within the time scale tq around the previous rate 
change than the time scale longer than tq. Such a case could 
be possible if the system provides some redundant information 
about the point of the rate change to customers as quickly as 
possible. 

For this choice of the observation time distribution Pq (t) — 
e~*^^°, we have 



<5„ = — 

"^0 Jo 



drPw (t) e 



(24) 



Therefore, from equations (fT2l l. (fl4l i and (fTSl l. it is possible for 
us to derive the waiting time distribution i7(s) and average 



waiting time and the deviation from the value. For the choice 
of a Weibull first passage time distribution Pw{t), we need 
some algebra, but easily obtain 

e«/-V^(To;s) 



However, at the same time, the standard deviation also behaves 
as 



(25) 



fli/" ^r(l/TO)r(3/m) - r(2/m) 



-O(V'). (31) 



where we defined ^^{tq; s) by the following integral form 



dzexp 



To 



(26) 



In Fig. [8] we plot the distribution for the case of several values 
of To and m — 0.585 and a — 49.6345. From this figure, we 
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Fig. 8. The waiting time distribution Q(s) for tlie case of the weight function 
^o(*) = c"'/'^" with several values of tq. For each plot, we set m = 0.585 
and a = 49.6345. 

find that the probability for large waiting time s increases as 
the relaxation time tq increases. We should notice that the 
case of To ^ oo turns out to be random observation by the 
customers. 

From the equations ( fT4l i. the first two moments are obtained 

by 

w =^ (s) — 



This means that even if the trader observe the rate according to 
a priori knowledge, namely, the observation time distribution 
Po{t) = c^*/^", the standard deviation from the average 
waiting time is the same order as the average waiting time 
itself. 

In this paper, we considered just only two cases of the 
observation time distribution Po{t), however, the choice of 
the distribution is completely arbitrary. The detail and more 
carefully analysis of this point will be reported in our forth 
coming article [17]. 

3) Comparison with empirical data analysis: It is time for 
us to compare the analytical result with that of the empiri- 
cal data analysis. For uniform observation time distribution 
Po{t) = 1> we obtained a — 60.2284 minutes. On the other 
hand, from empirical data analysis, we evaluate the quantity 
(O by sampling the moment as £:(t") = (1/A^)E^it" 
directly from Sony bank rate data [16] and find a = 74.3464 
minutes. There exists a finite gap between the theoretical pre- 
diction and the result by the empirical data analysis, however, 
the both are the same order. The gap might become small if 
we take into account the power-law tail of the first passage 
time distribution. In fact, our preliminary investigation shows 
that for the average waiting time, the power-law tail makes the 
loooogap between the theoretical prediction and empirical observa- 
tion small [17]. Therefore, the same tail-effect is reasonably 
expected even in the analysis of the deviation. 



where we defined 
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It should be noted that for large tq, the leading order of the 
function ^m,o(T) is evaluated as ^m,o(T) = 1 + 0{t^^). Thus, 
the average waiting time w for Po{t) = e"*/"^" becomes 
smaller than that obtained for Pq (t) = 1 as 

r(-) 



W 



A/rn_ 



r(-) 



(30) 



VI. Concluding remarks 

In this paper, we proposed a different procedure from 
the conventional derivation of the renewal-reward theorem. 
This derivation enables us to evaluate arbitrary order of the 
moments of the waiting time of the on-line foreign exchange 
trading rate for the individual customers. We directly derived 
the waiting time distribution and evaluated the deviation 
around the average waiting time, which is not supported by 
the renewal-reward theorem, for the Sony bank USD/JPY 
exchange rates. We tested our analysis for several cases of 
the observation time distribution of the customers and found 
that the average waiting time and deviation from the value 
are the same order even if the customers possesses a priori 
knowledge as a form of the observation time distribution as 
Po{t) = This result might be understood as follows. 

As we mentioned, the system we dealt with in this paper has 
two types of fluctuations, namely, fluctuation in the intervals of 
events (price changes) and fluctuation in the observation time 
for the individual customers to observe the rate through the 
World Wide Web. As well-known, in the A^-body systems {N 
is extremely large), there exists theremal fluctuation (quantum- 
mechanical fluctuation as well in low temperature) affecting on 
each particle, however, the macroscopic variables (the averages 
over the Gibbs distribution) like pressure or internal energy 



are determined as of order 1 objects and the deviation from 
the average becomes zero as This is a reason why 

statistical mechanics could predict a lot of physical quantities 
and we could compare the prediction with the same quantity 
which was observed in experiments. On the other hand, 
although the financial market price changes are taken place 
as a result of trading by many people, the observed rate itself 
is regarded as a result of effective single particle problems (a 
complicated single random walker) and there is no such scale 
parameter like the number of particle A^. This is an intuitive 
account of the result, namely, the average and the deviation 
are the same order even if there exist two kinds of fluctuation 
in the systems. 

Remaining problems concerning the analysis of the Sony 
bank rate are firstly to investigate the effect of the tail of the 
first passage time distribution. Our preliminary observation 
from the empirical data implied that, at some critical point, 
the distribution changes its shape from a Weibull-law to a 
power-law [26]. Therefore, it is important for us to check 
to what extent the prediction by the renewal-reward theorem 
is modified by taking into account the tail effect of the first 
passage time distribution. Secondly, we should show explicitly 
that the first passage time distribution described by the Mittag- 
Leffler function [9], [10], [15]: 



P(r) = 



n=0 



f3n 



jr/ro) 
r(/3n+l) 



(0 < /? < 1) (32) 



is impossible to realize the first passage process of the Sony 
bank by comparing the analytical result of the average waiting 
time or the Gini index with those obtained by the empirical 
data analysis. The detail of these studies will be reported 
shortly in [17]. 

As we showed in this paper, our queueing theoretical 
approach might be useful for us to build artificial markets 
such as the on-line trading service so as to have a suitable 
waiting time for the individual customers by controlling the 
width of the rate window. Moreover, theoretical framework we 
provided here could predict the average waiting time including 
the deviation. 

We hope that our strategy in order to analyze the stochastic 
process of markets from the view point of the waiting time of 
the customers might help researchers or engineers when they 
attempt to construct a suitable system for their customers. 
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